Frontiers in Molecular Biosciences

Background: Stress-induced childhood-onset neurodegeneration with variable ataxia and seizures (CONDSIAS) is a rare autosomal recessive disorder caused by biallelic variants in ADPRHL2, which encodes ADP-ribosylhydrolase 3 (ARH3), a key enzyme involved in poly (ADP-ribose) (PAR) metabolism. Although Minocycline has been reported to attenuate PAR-mediated neurotoxicity primarily through modulation of PARP-dependent pathways, whether it may also interact with ARH3 or influence the structural behavior of pathogenic ARH3 variants remains unknown. This study was designed to explore this possibility by integrating clinical observation with computational structural analyses. Methods: Comprehensive clinical evaluation, targeted Sanger sequencing, and in silico pathogenicity analyses were performed. Protein modeling, molecular docking, and 100-ns molecular dynamics simulations were conducted to evaluate the predicted structural consequences of the p.Thr79Pro variant and to explore potential interactions between ARH3 and Minocycline. Results: A homozygous ADPRHL2 variant (NM_017825.3:c.235A>C; p.Thr79Pro) was identified in a child with CONDSIAS. Computational analyses predicted reduced structural stability and increased conformational flexibility of the mutant ARH3 protein relative to the wild-type structure. MM-GBSA calculations estimated differences in binding free energies between the wild-type (-34.51 kcal/mol) and mutant (-39.76 kcal/mol) ARH3-Minocycline complexes, suggesting subtle differences in their predicted energetic profiles. Clinically, neurological progression appeared stable, with improved motor function observed during approximately one year of follow-up and no notable treatment-related adverse effects. Conclusions: By integrating clinical observations with computational structural analyses, this study provides preliminary computational support for the hypothesis that Minocycline may influence ARH3 conformational behavior in addition to its proposed effects on PARP-dependent pathways. Although these findings do not demonstrate direct molecular binding or therapeutic efficacy, they provide a biologically plausible framework for future biochemical, cellular, and functional investigations. Keywords: CONDSIAS; ADPRHL2; ARH3; Minocycline; molecular docking; molecular dynamics simulation; structural bioinformatics; translational medicine

The ability of small molecules to cross the blood-brain barrier (BBB) remains a major bottleneck in neurotherapeutic development. While experimental assays and machine learning approaches provide approximate permeability estimates, they lack atomistic insight into the underlying transport mechanisms. Here, we employ all-atom molecular dynamics simulations of a compositionally realistic BBB lipid bilayer to characterize the passive permeation of two bioactive propolis-derived compounds, Caffeic Acid Phenethyl Ester (CAPE) and Artepillin-C (ARC). Using steered molecular dynamics and umbrella sampling, we computed free energy profiles, diffusion coefficients, and permeability metrics across the membrane. CAPE encounters a modest barrier at the lipid headgroup region but minimal resistance within the hydrophobic core, resulting in a low free energy barrier ([~]2-3 kcal/mol) and favorable permeability (logP_eff {approx} 0.28). In contrast, ARC exhibits a substantial energetic barrier within the membrane core, leading to high resistivity and strongly unfavorable permeability (logP_eff {approx} -10.91). The heterogeneous lipid model reproduces experimentally consistent membrane properties and reveals how lipid composition modulates transport energetics. These findings provide mechanistic insight into BBB permeability and demonstrate the utility of atomistic simulations for guiding the design of neuroactive therapeutics.

Genome sequencing of the heterogeneous primary mitochondrial disorders (PMD) frequently reveals variants of uncertain significance that require functional tests for diagnosis, and does not identify variants in all patients. We analyzed mitochondrial enzyme assays, blue native polyacrylamide gel electrophoresis (BN-PAGE) with in-gel activity staining, complex I assembly blot, and select protein abundances in fibroblasts of a case series of 204 PMD patients divided into functional classes, in comparison to 51 controls and 53 differential diagnostic conditions. Overall, sensitivity and specificity for respiratory chain enzyme assays were 46% and 93% respectively, for BN-PAGE 40% and 98%, for complex I assembly assay 49% and 99%. The overall sensitivity of all tests was 76%, specificity 93%, with positive predictive value 96% and negative predictive value 67%. Categories with high sensitivity were isolated complex deficiencies, nuclear DNA-encoded mitochondrial protein synthesis defects, co-factor defects, and mitochondrial amino-acyl-tRNA synthetase conditions when aided by protein abundance. Mitochondrial DNA mutations and maintenance disorders showed poor sensitivities. Secondary dysfunctions were rare. A complete battery of functional tests showed strong diagnostic clinical utility in fibroblasts.

Emerin (EMD), an inner nuclear membrane protein essential for nuclear architecture integrity, gene expression, cellular signaling, and chromatin stability, interacts with the LINC complex and participates in cytoskeleton-nucleoskeleton communication by binding to nuclear actin filaments. EMD is implicated in migration, invasion, and metastasis in some tumors, but its role in glioblastoma (GBM) remains unclear. This study evaluated the effects of EMD knockdown and overexpression in GBM cell lines following genotoxic treatment with cisplatin. In both wild-type p53 (U87-MG) and mutant p53 (U138-MG) GBM cells, EMD expression is high, and cisplatin treatment did not affect these protein levels. EMD knockdown in U87-MG cells significantly increased cisplatin IC50, viability, and proliferation. Conversely, stable overexpression of EMD in U87-MG cells led to reduced cisplatin IC50, viability, proliferation, and migration. EMD knockdown or overexpression did not affect any U138-MG phenotypes, with or without cisplatin treatment. Modulation of EMD levels causes morphological changes in stress fiber cytoskeleton, whereas overexpression of EMD in U87-MG cells promotes an increase and a decrease in nuclear and cytoplasmic actin levels, respectively. These biological responses of U87-MG cells overexpressing EMD were coincidentally associated with alterations in the levels of pH2AX(Ser139), p-p53(Ser15), p53, and p21Kip1 proteins after cisplatin exposure. In sum, modulation of EMD levels affects the viability, migration, and proliferation of wild-type p53 GBM cells treated with cisplatin, suggesting unknown roles in the DNA damage response and repair. This work highlights EMD as a potential regulator of GBM chemoresistance and a target for therapeutic intervention.

Protein behavior inside cells is dominated by the crowded nature of the intracellular environment. Progress in structure determination of proteins and protein complexes, based on advances in Artificial Intelligence, provides an opportunity for structure-based modeling of cellular phenomena. Such modeling at the atomic resolution has been advanced by the traditional simulation techniques, e.g. molecular dynamics. A recently developed docking-based approach implements Markov Chain Monte Carlo sampling of intermolecular energy landscapes, offering several orders of magnitude faster simulation protocols. The approach allows addressing much longer trajectories of macromolecular systems in the crowded intracellular environment at atomic resolution. The sampling by design avoids low-probability (high-energy) states, which greatly accelerates the simulation process. A notable feature of this docking-based approach is the rigid body approximation of protein structures. The rigid-body approximation had been the primary direction in the protein docking field up until recent developments in deep learning. The rigid-body approach should be quite robust for the higher energy transient interactions that dominate the highly crowded cellular environment, as they likely involve relatively small conformational change. However, it is less applicable to the low-energy protein-protein complexes, especially those involving flexible regions. We addressed this problem by incorporating AlphaFold3 top models of the protein complexes in the mapping of the intermolecular energy landscape, as representative of the low-energy configurations of the protein assembly. By the nature of the AlphaFold predictions, these models involve appropriate conformational change between unbound and bound structures. These low-energy docking poses are combined with the rigid-body docking predictions that cover the multiplicity of the transient interactions. Such combination directly addresses the conformational flexibility of proteins upon binding along with the multiplicity of the transient protein encounters in the crowded cellular environment. SIGNIFICANCEProtein behavior inside cells is dominated by the crowded nature of intracellular environment. A recently developed approach allowed addressing long simulation trajectories of macromolecular systems in such environment at atomic resolution. A notable feature of this approach is the rigid body approximation in representation of the protein structures, which had been popular in the field up until the recent developments in artificial intelligence. However, such approximation is less applicable to stable protein-protein complexes, especially those involving flexible regions. We addressed this problem head-on by incorporating top deep learning-generated models of protein complexes. The new approach directly accounts for the flexibility of protein structures upon binding, along with the multiplicity of the transient protein encounters in the crowded cellular environment.

A priori assessment of target proteins' druggability remains an unsolved problem in the field of drug development. The empirical approaches widely used to solve this problem demonstrate low efficiency. In this work, we investigated the factor of hydration of a representative set of 65 evolutionarily and structurally unrelated human enzymes in a water environment. This factor depends only on the structure of the proteins, and not on the physical and chemical properties of any potential ligands. The results show that, unlike the widely used approaches based on calculations of the accessible surface area (ASA), the content of low-entropy water molecules (LEW) in the active sites of human enzymes is systematically higher than that in other areas of their surface, including inactive cavities. Optimal criteria and a step-by-step procedure for identifying protein ligand binding sites are proposed. The proposed approach, based on the calculation of the LEW content in the first hydration layer of potentially interesting target proteins, makes it possible to evaluate their medicinal suitability even before the development of any ligands. The article also presents the results of a comparative analysis of experimental Raman spectroscopy data and the results of molecular dynamics simulations of water hydrogen bonds using three widely used water models (TIP3P, OPC3, and TIP5P) and standard algorithms for calculating hydrogen bond networks.

Significance Missense mutations can lead to pathological effects in human cells. Predictive methods that account for structural context, such as AlphaMissense, can provide pathogenicity scores. The accumulation of pathogenicity hotspots can reveal important structural features within individual proteins of protein families, such as GLUT transporters. Mapping pathogenicity scores onto the structure can thus provide a mechanistic explanation of the protein function necessary for its role in the cell. Abstract Non-synonymous amino acid substitutions (missense mutations) are common in the general population; some are causative of serious disease. Depending on their structural context, they can disrupt protein function, folding, or dynamics. Computational predictive methods developed in recent years, such as AlphaMissense, provide new insights into how missense mutations affect protein structure by predicting and mapping their pathogenicity across each amino acid in the human proteome. In this study, we identify recurring patterns of pathogenicity prediction across the GLUT family membrane transporters encoded by genes slc2a1-14. Within the GLUT transporter family, we observe higher pathogenicity profiles in the transmembrane domains, particularly in pore-lining and binding-site residues. Predicted missense pathogenicity is elevated throughout residues assigned to the central cavity, suggesting sensitivity of the transport pathway. Another finding shows higher pathogenicity in specific transmembrane helices of the protein, with the same pattern across all proteins. On the other hand, we observed lower pathogenicity values in some representatives of the GLUT family. These findings show that the pathogenicity of glucose transport within the GLUT family may be shaped by functional redundancy and physiological essentiality across GLUT groups.

DnaK, a prokaryotic Hsp70 chaperone, plays a central role in proteostasis by restoring native structures to heat-denatured proteins in an ATP-hydrolysis-dependent manner. While structures of DnaK in complex with nucleotides, co-chaperones, and short peptides have been resolved, structures with larger, stably folded substrates--such as firefly luciferase (Fluc, 61 kDa)--are lacking, limiting mechanistic understanding of how DnaK refolds such proteins. Here, we generated models of the DnaK-Fluc complex using AlphaFold3 and evaluated their mechanistic relevance. In one of three major model clusters, Fluc is unexpectedly immobilized beneath the DnaK -helical lid against the nucleotide-binding domain (NBD), rather than interacting primarily with the substrate-binding domain {beta} (SBD{beta}), as commonly assumed. All-atom molecular dynamics simulations indicate that, in this configuration, the lid can engage a thermally destabilized Fluc helix (residues 405-411), which we recently identified as the first--and likely the only--helix to irreversibly melt at 42 {degrees}C. Upon binding, the lid forms extensive hydrogen-bonding interactions with the melted helix. These interactions persist during lid movement toward SBD{beta} (following ATP hydrolysis), enabling the lid to actively extract the helix from the Fluc surface. In contrast, simulations with the helix in its native folded state show that the lid cannot extract it, leaving the native structure unaffected. Equilibrium simulations further indicate that, once extracted and mechanically stretched, the melted helix can refold to its native conformation. Together, these findings suggest a revised mechanism for DnaK-mediated protein refolding, in which the -helical lid selectively recognizes structurally compromised segments, forms stabilizing hydrogen bonds, and--powered by ATP hydrolysis--mechanically pulls them away from the protein surface to facilitate their refolding. SIGNIFICANCEDnaK is a model chaperone, which can reactivate thermally denatured proteins. Over the span of 40 years, significant findings have been made about DnaKs structure, dynamics and interactions with its co-chaperones, the exact molecular mechanism by which DnaK refolds misfolded proteins remains a mystery. This work exploited Alphafold3 to generate atomistic models of complexes between DnaK and Firefly luciferase. Molecular dynamics simulations directly captured how DnaK may assist thermally denatured proteins by mechanically pulling out their misfolded helices. This study provides a new insight into the DnaK mechanism.

Bordetella pertussis, the causative agent of whooping cough, is a reemerging public health threat. While the Tripartite Tricarboxylate Transporter (TTT) system BctCBA was previously implicated solely in citrate uptake, we demonstrate that the solute-binding protein BctC specifically binds citrate chelated with Zn{superscript 2} and Ni{superscript 2}. To elucidate the molecular mechanism of this interaction, we determined the crystal structures of BctC in three states: apo, open, and closed (citrate-zinc-bound), defining the structural determinants for metal-citrate recognition. Comparative analyses suggest that citrate-mediated divalent cation binding is a widespread feature among bacterial TTT homologs. Finally, in silico modeling of the full BctCBA complex predicts an elevator-type transport mechanism. Together, these findings redefine the functional scope of BctCBA, revealing a sophisticated strategy by which B. pertussis exploits organic chelators to acquire essential trace metals during infection.

Prestin (SLC26A5), a membrane protein in cochlear outer hair cells, drives electromechanical transduction essential for mammalian hearing. Unlike other SLC26 anion transporters, prestin functions as a voltage-dependent molecular motor, transitioning between compact and expanded conformations. How this transition relates to the transporter cycle of SLC26 family members remains unclear. Here, multi-microsecond molecular dynamics simulations starting from the compact state reveal a rapid, spontaneous transition to an expanded state that resembles the inward-facing conformation of the anion exchanger pendrin (SLC26A4 from mouse). An accompanying transmembrane area expansion is localized to the inner membrane leaflet, likely leading to membrane bending. In line with this observation, reduced unitary sensor charge movement accompanies neutralization of charged residues localized near the inner leaflet. Simulations also uncover a previously uncharacterized compact conformation resembling outward-facing pendrin and predict an extracellular anion-binding site in prestin. In fact, in the presence of thiocyanate anions, we observe a previously unresolved binding site in a 3.27-[A] cryo-electron microscopy structure of prestin. Furthermore, like prestin, pendrin exhibits a non-linear capacitance, an indication of voltage-dependent conformational switching. Together, these findings suggest that prestin and pendrin share core structural and functional properties, notably parallels between expansion-contraction states and transporter function, though transition speeds may differ.

BackgroundMetabolic reprogramming is a common occurrence in tumor cells, where enhanced glycolysis promotes cell growth, invasion and migration. NF1 is tumor suppressor gene that downregulates the encoded neurofibromin protein. However, the effects of NF1 on energy metabolism and epithelial-mesenchymal transition (EMT) in glioblastoma multiforme (GBM), as well as the underlying molecular mechanisms, remain unclear. MethodsCRISPR/Cas9 gene editing technology was employed to construct GBM cell lines with NF1 gene mutations. Metabolomics was utilized to examine the impact of NF1 on metabolic remodeling in GBM. The Seahorse XF24 extracellular flux analyzer was used to detect the effect of NF1 knockdown on glycolysis and mitochondrial oxidative phosphorylation in GBM cells. Wound healing assay and Transwell chamber assay were utilized to detect the effect of NF1 on GBM cell invasion. Orthotopic tumor model in nude mice was established to explore the role of NF1 in vivo. In addition, Co-IP, western blotting, and immunofluorescence were used to explore the changes of key enzymes in glycolysis and mitochondrial oxidative phosphorylation and the relationship between NF1 and MFN1. ResultsThe expression of NF1 is decreased in glioma tissues and is significantly correlated with patient prognosis. NF1 knockdown may promote the invasion, migration, and EMT of GBM cells. At the same time, the activation of the AKT/mTOR signaling pathway promotes aerobic glycolysis in GBM cells, promotes mitochondrial division through targeted regulation of MFN1, and inhibits mitochondrial oxidative phosphorylation. NF1 deficiency promotes EMT in GBM cells by enhancing aerobic glycolysis and mitochondrial division. ConclusionNF1 deficiency promotes GBM glycolysis by activating the AKT/mTOR signaling pathway and inhibits the mitochondrial oxidative phosphorylation by regulating MFN1; NF1 deletion promotes GBM EMT by remodeling the pattern of energy metabolism.

Molecular dynamics simulations of nucleic acids are performed using a solvent-buffer distance of 10 [A] between the solute surface and the simulation box boundary. Although this cell size has been extensively explored in protein simulations, its implications for nucleic acid dynamics are not well understood. Nucleic acids are elongated, highly charged, and flexible structures with hydration and dynamical properties distinct from those of proteins and therefore, they may require different solvent-layer considerations in simulations. In this study, we investigated the effect of simulation cell size on nucleic acid dynamics by simulating a 30-base-pair double-helical nucleic acid structure and its two single-stranded forms using solvent-buffer distances of 3, 5, 10, 15, and 20 [A]. Smaller cells may impose restricted hydration, molecular crowding, and periodic image interactions. However, larger cells provide solvent space for conformational relaxation. A total of 45 s of molecular dynamics simulations were performed (3 structures x 5 cell sizes x 3 replicates x 1 s). Our results show that while the commonly used 10 [A] buffer may be sufficient to maintain the stability of the double-stranded nucleic acid, larger cells are required to capture the conformational dynamics of single-stranded structures. In both, increasing the cell size to 15 or 20 [A] enables broader conformational sampling. The first hydration shell exhibits reduced crowding in the 20 [A] cell, consistent with more relaxed conformations. At larger cell sizes, single-stranded nucleic acids adopt compact, self-associated conformations for stability. Together, this study presents physical insight into how simulation cell size and solvent environment influence nucleic acid dynamics.

Heteromeric amino acid transporters (HATs) mediate essential amino acid flux across membranes, but the molecular dynamics of substrate translocation remain poorly defined for many family members. Here, using conventional and adaptive steered molecular dynamics (cMD and ASMD) simulations, we identify residue W230 in the b0,+AT transport channel as a dynamic gate that regulates arginine (Arg) influx through side chain flipping. By integrating dynamic network analysis with dynamical cross-correlation of residue motions, we show that regulatory signals propagate from the Arg binding site through transmembrane helix 5 (TM5), a connecting loop, and TM6 to reach W230. We propose a dynamic gating mechanism for b0,+AT - mediated amino acid transport. Arg binding at V186 triggers signal propagation that enhances cooperative interactions between W230 and Arg, driving the side chain flipping of W230. Our findings reveal a dynamic gating mechanism underlying b0,+AT - dependent Arg transport and suggest that residue-triggered side chain reorientation may represent a conserved and efficient strategy in transporter function. Author SummaryAmino acids are the essential building blocks of life, and their transport across cell membranes is vital for nutrition and cellular signaling. Heteromeric amino acid transporters (HATs) mediate this process, yet how they physically move substrates through the protein at the atomic level remains poorly understood. In this study, we used advanced computer simulations to observe, in unprecedented detail, how b0,+AT--a key HAT member--transports the amino acid arginine. Our simulations revealed that a single residue, tryptophan 230 (W230), functions as a molecular gate: its side chain flips open to allow arginine to pass and then closes behind it, ensuring one-way traffic into the cell. We further discovered that the initial binding of arginine sends a signal through specific structural elements (helices and loops) to trigger this gate opening. This work not only uncovers a dynamic gating mechanism for b0,+AT but also suggests that similar side-chain flipping events may represent a common and efficient strategy used by other transporters to control substrate movement. Our findings provide a new framework for understanding transporter function and could inform future drug design targeting these critical membrane proteins.

The eukaryotic oligosaccharyltransferase (OST) is the enzyme responsible for initiating N-glycosylation of secreted proteins by transferring a pre-assembled lipid-linked oligosaccharide (LLO) donor to target asparagine residues most often found within N-x-S/T consensus sequences, or sequons. OST preferentially selects LLO donors with a distinctive glucoside Glc-(1-2)-Glc-(1-3)-Glc-(1-3)-capping the A-branch. After the N-glycosylation reaction, this motif is cleaved in a stepwise manner from the immature N-glycan structure before the folded glycoprotein exits the endoplasmic reticulum quality control (ERQC) cycle. While the -Glc-(1-3)-Glc-(1-3)-disaccharide is an important flag regulating binding to the calreticulin/calnexin chaperones, the terminal Glc-(1-2)-is removed immediately after OST catalysis, suggesting that its biological function may be directly linked to the OST catalytic efficiency. To understand how and why this capping motif affects the OST N-glycosylation efficiency, we rebuilt 3D models of the yeast OST in complex with an acceptor peptide and LLO donors substrates with and without terminal Glc-(1-2)-, and analysed their stability and dynamics with all-atom molecular dynamics (MD) simulations through both conventional, and Gaussian-accelerated (GaMD) sampling schemes. Our results indicate that the terminal Glc-(1-2)-is essential to anchor the full-length LLO donor to the OST through a complex network of intermolecular contacts extending from the catalytic site to distal subdomains. We show how this contact network is crucial to preserve the LLO catalytically productive alignment of its reducing end. We also show that the removal of the terminal Glc-(1-2)-leads to an increased flexibility of the LLO, which displaces the reducing end and redistributes the conformational ensemble towards misaligned states, which are less catalytically productive. These results provide a mechanistic basis linking the catalytic efficiency of the eukaryotic OST to the distinctive glucosylated structure of the LLO donor.

Mechanistic elucidation and inhibition of the pathogenic aberrant mitochondrial localization of the RNA/DNA-binding protein, TDP-43, can help in the therapeutics of the neurodegenerative disease amyotrophic lateral sclerosis (ALS). A mitochondrial localization sequence of TDP-43, M1, is largely solvent inaccessible, therefore, how it interacts with the mitochondrial import machinery to facilitate TDP-43s transit to mitochondria is unclear. Towards this, we examined the unfolding TDP-43s N-terminal domain (NTD) that hosts M1, using equilibrium all-atom molecular dynamics (MD) simulations, and observed an early loss of the hydrogen-bonded interactions between {beta}4-{beta}5 bridge and the interactions involving residues Phe-35 and Gly-40 of M1, indicating structural lability of M1 to become solvent-accessible that may enhance its interaction with the mitochondrial receptor(s) for import. Furthermore, via virtual screening of 2,115 FDA-approved and 515,545 non-FDA-approved small molecules from ZINC15 database towards binding to M1 and inhibiting TDP-43s mitochondrial import, we identified a molecule, ZINC73240059, that was previously characterized as an inhibitor of MAP kinase-activating protein kinase 2 (MAPKAPK2). ZINC73240059 remains stably bound to M1 of NTD during MD simulations manifesting negative Gibbs free energy ({Delta}G) with significant contribution from Pro-36 of M1. Overall, ZINC73240059 can be a molecule of interest towards thwarting TDP-43s pathogenic mitochondrial localization in ALS.

Malignant brain tumors are among the most aggressive and difficult to treat human cancers. Glioblastomas (World Health Organization grade IV gliomas) are particularly lethal and refractory to treatment. Few drugs exist that are even somewhat effective. Our investigation of the physiologic role of fatty acid (FA) activating enzymes (acyl-CoA synthetase; ACS) identified an ACS that was widely expressed in gliomas but not in normal glial cells. Depletion of this enzyme, ACSVL3 (very long-chain ACS3), by knockdown or knockout decreased the malignant behavior of several glioma cell models including U87MG and Mayo-22 cells both in culture and when grown as xenografts. Hypothesizing that ACSVL3 is a potential therapeutic target in glioma, we conducted a search for inhibitors of this enzyme and found that CB5 (grassofermata) was a promising candidate. Treating U87MG glioma cells with CB5 slowed growth in monolayer culture; the growth rate was similar to that seen in cells in which ACSVL3 was either knocked down or knocked out. CB5 inhibited growth in a dose-dependent manner over a narrow range, and concentrations above 10 M were toxic. Treatment at the lower dose of 3 M inhibited growth of U87MG cells but was reversible, suggesting that this dose was not toxic. CB5- treated U87MG cells exhibited an altered morphology with a larger size and longer projections. In contrast, normal human fibroblasts treated with 10 M CB5, a concentration that was toxic to U87MG cells, showed no effect on either growth rate or morphology. Treating U87MG cells with 3 M CB5 induced differentiation as shown by increased expression of the astrocyte-specific marker glial fibrillary acidic protein (GFAP). In contrast, GFAP levels remained low in ACSVL3 knockdown cells. CB5- treated U87MG cells were less invasive, and thus less malignant, than either untreated cells or ACSVL3 knockout cells when assessed by a scratch wound healing assay. Acute treatment of U87MG cells with 3 M CB5 decreased the ability of these cells to degrade FA of differing chain lengths from 16-24 carbons by {beta}-oxidation, suggesting that decreased ACS enzyme activity contributes at least in part to the drugs mechanism of action. NOD/SCID mice receiving up to 32 mg/kg/day CB5 by intraperitoneal injection showed no obvious side effects, suggesting that the drug was well-tolerated. Xenografts induced by subcutaneous injection of U87MG cells in the flanks of NOD/SCID mice were allowed to grow for 8 days after which half of the mice were treated with 2 mg/kg/day CB5. After 7 days of treatment, xenograft growth slowed in the treated mice and by 12 days tumor size had begun to decrease, suggesting therapeutic efficacy. When a similar study was done using xenografts induced by subcutaneous injection of Mayo-22 cells, which are maintained as subcutaneous tumors in mice rather than in cell culture, the effect of CB5 on tumor growth or weight at sacrifice was not statistically significant. The results of these studies suggest that CB5 may have therapeutic value in malignant glioma. Additional studies using other glioma models and other drugs chemically related to CB5 seem warranted.

The limited understanding of transporter substrate spectra constrains our ability to interpret cell and membrane function, highlighting the need for methods that enable transporter deorphanization and characterization of promiscuous transport activities. Here, we present a Xenopus oocyte-based platform for unbiased transporter substrate discovery. Oocytes expressing heterologous solute carrier proteins (SLC) were incubated in human blood serum, a chemically complex metabolite library containing thousands of endogenous metabolites and xenobiotics, followed by paired untargeted LC-MS/MS profiling of intracellular extracts and surrounding medium to capture metabolite exchange events. Across the five human SLC transporters, viz. SLC10A2, SLC10A6, SLC13A2, SLC16A10, and SLC46A1, metabolite exchange signatures were detected, and automated feature annotation was refined by manual chromatographic peak inspection. The workflow recovered known substrates of SLC10A2 and SLC16A10 and identified additional transported metabolites with MS/MS confirmation. This method provides a scalable framework for transporter substrate profiling and prioritization of candidates for targeted validation.

Gain-of-function (GoF) variants in the Kv7.2 channel are associated with a clinically relevant subset of neurodevelopmental disorders. While most GoF substitutions identified so far affect the voltage-sensing domain, we recently described three mutations in the intracellular-facing activation gate (AG), G313S, A317T, and L318V. Electrophysiological recordings showed that these variants increase macroscopic current density and enhance channel open probability. Consistently, molecular dynamics (MD) simulations revealed that they hinder complete channel closure by destabilizing the closed AG and increasing hydration of the central cavity (CC). Whether this partially open conformation can support K+ permeation, however, remained an open question. Here, we combined long-timescale atomistic simulations and simulations with applied electric fields to evaluate the stability of these mutant-associated AG states over longer timescales and their functional relevance. In new trajectories, all three substitutions consistently shifted the closed intracellular gate toward a widened, water-accessible conformation, accompanied by increased CC hydration. We then assessed the functional significance of this partially open state by simulating the A317T channel under applied electric fields. The conformation supported K+ translocation, whereas the closed WT pore remained impermeable under all tested voltages. When simulations were started from open channel conformations, both WT and A317T conducted K+ ions, indicating that the main effect of the substitution is to destabilize closure of the intracellular gate rather than to alter the fully conductive open state. Together, these data show that AG GoF variants can generate an intermediate gate conformation that permits ion permeation, providing a mechanistic link between mutant-induced pore remodeling and Kv7.2 dysfunction in KCNQ2-related disease.

Magnetic treatment in tomato seeds (Solanum lycopersicum L.) has been studied as a biotechnological technique to induce a reduction in germination times and enhance plant development. However, the modified cellular mechanisms involved in the reduction of germination times or the improvement of development parameters are not yet clearly established. To explore a possible altered cellular mechanism, the effect of homogeneous static magnetic fields on the structure of the cyclic nucleotide-gated channel 6 (CNGC6), the modification in the organization of POPC lipids in the plasma membrane, and changes in calcium ion mobility were evaluated. For this purpose, coarse-grained molecular dynamics simulations were performed using the Martini 3 model in GROMACS, applying five different magnetic flux densities (0.000, 0.001, 0.010, 0.100, 1.000, and 10.000) T over 1 000 ns. The results showed an anisotropic effect in the longitudinal direction of the protein, which generated heterogeneous behavior among the chains of the homotetramer; this altered the conformation of the CNGC6 channel and modified the pore bottleneck. In contrast, no significant changes were observed in the conformational order of the POPC phospholipid chains. As a preliminary, single-replicate exploratory study, these results suggest that homogeneous static magnetic fields may induce specific structural modifications in the CNGC6 ion channel of Solanum lycopersicum L. without compromising the integrity of the lipid bilayer or the dynamics of ion transport within the analyzed timescale; these preliminary findings provide a molecular-level structural basis for future experimental and computational investigations of magnetic field effects on plant cyclic nucleotide-gated channels.

Amylose contributes to starch crystallinity, but the stability of packed amylose double helices in water at elevated temperature remains insufficiently characterized. Here, we used molecular dynamics simulations to test whether chain length affects the short-timescale stability of A-type amylose oligomers in water. Six systems differing in chain length (6, 12, or 24 glucose units per chain) and oligomer size (isolated double strand or dodecamer of six double strands) were simulated, and five independent 1-s production runs were analyzed for each simulated condition. Oligomers with six glucose units showed structural collapse accompanied by increased water penetration. By contrast, dodecamers with 12 or 24 glucose units largely retained packed double-helical organization over the simulated timescale, although fraying was observed at their ends. These results indicate that chain length and lateral packing strongly affect the early structural response of amylose-like crystalline segments in hot water. The present simulations do not establish the ultimate fate of longer oligomers at longer timescales, but they identify a relative stability difference that is relevant to molecular interpretations of hydration-driven disordering in starch.